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Abstract 

We study the response of a recently proposed global coupling model of charge- 
density waves to a joint ac+dc drive. The model is the standard Fukuyama- 
Lee-Rice model with an additional global coupling term to account for in- 
teraction with (uncondensed) quasi-particles. We find that traditional mode- 
locking is accompanied by hysteresis on mode-locked steps. The experimen- 
tally observed asymmetry in the size of hysteresis is reproduced and explained 
by drawing attention to the importance of the pinning force. For large values 
of global coupling, increased rigidity leads to vanishing of subharmonic steps 

consistent with experiments. 
PACS numbers: 72.15.Nj, 05.45. +b 
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Charge- density wave (CDW) conductors have been established as ideal systems in which 
to study the dynamics of many degrees of freedom in the presence of disorder CDWs 
display a variety of nonlinear nonequilibrium phenomena, including dynamical phase transi- 
tions and mode- locking ||. In the presence of an applied dc- field Ed c , above a threshold 
E t h, the CDW depins and carries an average current J c d w which varies continuously with 
Edc- A similar behavior is observed at the mode-locking transitions, where in the presence 
of joint ac+dc drive the "narrow-band noise" frequency, u n b n oc J c dw, which characterizes 
the CDW current oscillations, locks onto the external frequency u ex with uj n b n /oj ex = p/q, 
p and q integers. The mode-locked state is the analog of the pinned state, with the CDW 
remaining locked at the same frequency for a finite range of E^. This defines a mode-locked 
step over which the average CDW current remains fixed. Upon unlocking at lower and 
higher "thresholds" the current varies continuously with Ed c . 

Over a range of temperatures the depinning and mode-locking transitions are smooth 
and can be analyzed in terms of dynamical critical phenomena 0. However, at lower 
temperatures, CDWs also display discontinuous and sometimes hysteretic transitions. In 
the case of the depinning transition, hysteresis, or "switching," is observed at a second, 
higher threshold field E* h More recently, Higgins et al. []5] observed hysteresis on mode- 
locked steps. Most notably, the size of the hysteresis loops on opposite edges of a given step 
exhibit a marked asymmetry: for certain values of the ac-drive the right edge of the step 
is hysteretic while the left edge displays a smooth mode-locking transition. The observed 
hysteresis is small in size and occurs in a parameter regime in which the sample is not 
switching in the absence of ac-drive. In a different set of experiments, Levy et al. || pointed 
to the chaotic behavior on the mode-locked steps in switching samples, in which case no well- 
defined mode-locking hysteresis loop can be identified. Instead, a period-doubling route to 
chaos with an associated low- dimensional attractor was observed while at lower frequencies 
high dimensional and noisy dynamics were seen. 

These hysteretic and chaotic phenomena provide new challenges to our theoretical under- 
standing of these systems. In particular, motivated by these experiments, we have studied 
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the mode-locking behavior of a new global coupling model Jj]] of CDW dynamics which 
has been seen previously to show some promise in explaining the temperature dependent 
behavior of CDWs by producing hysteresis at the depinning transition as well as delayed 
conduction ||. In this paper we focus our attention on the one- dimensional version of this 
model, which we find produces some of the experimentally observed results and also allows 
us to make further predictions; however, the Id model comes up short on some important 
aspects of the experiments as well. 

Many aspects of CDW dynamics including the phenomena of mode-locking are well 
described by the Fukuyama- Lee- Rice (FLR) model |J . An important property of this model 



is the so-called no passing rule |IJ which requires that the current be a single- valued function 
of Ed c . Hence the FLR model cannot describe the low-temperature hysteresis exhibited at 
the depinning or mode-locking transitions. As pointed out by Littlewood |11 and recently 



shown by Levy et al. 0, the role of normal carriers excited across the CDW gap (Aqdw) 
becomes increasingly important at low temperatures. This adds a new global coupling term 
to the FLR equations of motion, altering the I-V characteristics significantly. 

The global coupling model is obtained by incorporating current conservation within the 
FLR equation of motion. Ignoring the displacement current, the conservation equation is 
readily integrated in Id |7j leading to the following equation suitable for numerical simula- 
tions, 

(7o + Ji)fa = kA 2 ( f) i - sm(2n(fa + &)) + E(t) 

'Yi N ■ 
iv i=l 

Here fa represents the CDW phase at the ith impurity site, 70 is the time constant for CDW 
relaxation (we set 70 = 1), and 71 is the normal carrier resistance, assumed to have an 
exponential temperature-dependence: 71 cx exp(A C Dw/kBT). The FLR pinning parameter 
is k, so that the two limits 2nk <C 1 and 2irk ^> 1 correspond to the strong and weak 
pinning limits, respectively. The f3i G [0, 1) are random variables due to disorder, and A 2 is 
the lattice Laplacian. The average CDW current is given by J c d w =< fa >x,t— ^<fc + Fpin 



where the second term is the time and space averaged pinning force 

F pin = - < sin(27r(4 + A)) >x,t ■ (2) 

We simulate Eq. (1) with E(t) = Ed c + E ac cos(u ex t) using periodic boundary conditions in 
one dimension. Fifth-order Cash-Karp Runge-Kutta |T3| is used to integrate the equations 



for iV varying from 10 to 1024 with various values of k. 

When 71 = the standard FLR model is recovered. For non-zero 71 the time constant for 
4>i is adjusted and the global coupling term is relevant. It was shown that the additional 
term gives rise to hysteresis in the regular (E ac = 0) I-V curve; however, the size of hysteresis 
can be very small for small 71, with the difference between the pinning and depinning fields 
scaling roughly as 7^ for 71 <C 1 [ 12 1 . For the mode- locking transition, we define the size 



of hysteresis Sh ys t to be the magnitude of the difference between the locking and unlocking 
fields, and examine variations in Sh ys t in different regimes. 

Figure 1 shows mode-locking (E ac 7^ 0) I-V curves for different values of 71 in the 
weak pinning regime (k = 1.0). For 7! = 0.0 harmonic (q = 1) and many high order 
subharmonic (g ^ 1) steps are produced. As 71 increases hysteresis emerges asymmetrically 
on harmonic and some subharmonic steps. With increasing ji subharmonic steps become 
smaller as harmonic steps grow wider. The left edges of the steps become rounded as the 
right edges become sharp and hysteretic. No chaotic behavior is observed. Note that within 
a hysteresis loop (e.g. right edge of the 0/1 step in Fig. 1(b)) two distinct mode- locked 
states, i.e. different p/q ratios, sometimes exist for identical set of parameters, one on the 
upper branch and one on the lower branch. This means that nondegenerate mode-locked 
attractors exist for this model, where the FLR model only exhibits degenerate mode-locked 
attractors displaying the same velocity ||10| . We emphasize that in order to obtain this 
much structure in the I-V curves it was necessary to go to relatively large values of k. In 
comparison, while the strong pinning case of k = 0.01 also displays asymmetric hysteresis 
on mode-locked steps, the I-V structure was much simpler, displaying harmonic and only 
p/2 subharmonic steps even for large system sizes. 



In analogy with the behavior for the depinning transition |7|], we expect the Sh ys t to 
increase with increasing 71. This is always observed for strong pinning, and in general this 
characterizes the basic trends. However, in the weak pinning case, there are sometimes 
exceptions occurring in cases where we observe subharmonic steps on the upper branch of 
the hysteresis loops. These additional steps within a hysteresis loop appear to increase Sh ys t 
as is seen by comparing Fig. 1(b) and 1(c). This feature disappears for small k or large 71. 

In general, the disappearance of subharmonic steps throughout the I-V curve with in- 
creasing 71 is consistent with experiments ||. In the model internal deformations of the 
CDW must be accompanied by the screening currents of quasi-particles in order to conserve 
charge locally. With decreasing temperature quasi-particle conductivity (I/71) decreases 
exponentially, inhibiting the "backflow" current. This freezes out the internal CDW defor- 
mations and results in a rigid CDW. It is well known that subharmonic mode-locking is a 
feature of many degrees-of-freedom dynamics and that this behavior disappears in single 
degree-of-freedom models ||. Here, the disappearance is due to low dimensional behavior 
in a high dimensional model, arising because of the dominance of the global coupling. 

Figure 2 shows the "phase diagram" in the weak pinning (k = 1.0) regime for 71 = 2.5. 
Figure 2 differs from the standard phase diagrams where the step width is plotted in the 
E ac vs. Edc plane. Here, in order to illustrate the existence and variations of hysteresis, we 
plot the full I-V curves for different values of E ac . For each mode-locked state the width 
of the step oscillates as E ac is varied. Within the FLR model, in the high velocity limit, 
the oscillations in harmonic mode- locked steps are described by Bessel functions ||14j| . In 



Fig. 2 the dashed lines connecting the unlocking transition points appear to be at least 
qualitatively consistent with the FLR case though we have not investigated the functional 
form in detail. The hysteresis size asymmetry is most pronounced in the first half of the first 
oscillations (low E ac ) where large hysteresis is observed. For larger E ac , as in the second 
oscillation of the 1/1 step, hysteresis becomes much smaller and symmetric on the opposite 
edges of the steps. The asymmetry occurs even for high velocity mode- locked steps, e.g. 4/1 
step (not shown here), at low E ac values. Therefore, asymmetry appears to be a low E ac 



but not necessarily a low velocity phenomenon. 

We suggest that an important variable governing the size of hysteresis within a given 
phase diagram is the pinning force at the edges of the mode-locked step. In Fig. 3 we show 
a schematic asymmetric mode-locked step labeling important quantities. By drawing the 
J c dw = E dc line through the step, one can visualize the magnitude of the average pining 
force, l-Fpinl, as the height difference between the two curves. One might naively expect that 
in the middle of the step the pinning force would be zero so that at the edges of the step 
the pinning force would take on some extremum values as allowed by the value of E ac . This 
is indeed the case for large E ac ; the second oscillations, for example, exhibit steps for which 
the value of the pinning force at the center of the step is close to zero, and the value of 
Fpfy £ is equal in magnitude, with opposite sign, at the two edges. However, within the first 
oscillation, the center of the steps does not occur at E dc = i.e. where F pin = 0. The 
centers of these steps lie at higher values of E dc (nonzero pinning force) and this ultimately 
gives rise to the hysteresis asymmetry of the steps. The center of these steps occurs at higher 
values of E dc since for small E ac the pinning force cannot change appreciably from its value 
when E ac = 0. To increase J cdw to the mode-locked value, J^, one must go to higher E dc . 
Therefore, for low E ac the steps are "shifted" to the right as is apparent from Fig. 2. 

The above argument suggests that the 0/1 step should always be very nearly symmetric, 
since for the 0/1 step the pinning force is zero at E dc = 0, which is always close to the center 
of the step ||15|| . Higgins et al. [||] did not look at the left edge of the 0/1 step so this remains 
to be verified. However, they did discover that for steps other than 0/1, hysteresis occurred 
for those values of E ac where the steps were the widest, i.e. middle of the first oscillation. 
This is consistent with Fig. 2 where the largest hysteresis occurs in the middle of the first 
oscillation. More generally, the wider steps tend to have a larger pinning force at the edges, 
and the larger the pinning force the larger the size of hysteresis. 

The observed qualitative trends in the relationship between the size of the hysteresis 
(e.g. Shyst = E T , 2 — Et 3 for the right edge hysteresis in Fig. 3) and the magnitude of the 
pinning force (Eq. (2)) suggests that it may be worth searching for some quantitative scaling 
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relationship. Restricting our attention to a particular phase diagram in the strong pinning 
regime (k = 0.01 and 71 = 1.0, and w ex = 2n), we have attempted to find such a relationship 
with limited success. In particular, a scaling fit between Sh ys t and F^ e appears to hold 
well for the symmetric steps. However, an extrapolation of this scaling breaks down in the 
low E ac limit where the asymmetry is most pronounced. In that case on the right edges of 
the steps, Sh ys t is even larger than what the extrapolation predicts, while Sh ys t is smaller 
(often zero) on the left edges. This result is indicative of a crossover occurring between the 
high and low E ac limits, which can be accounted for in a crude manner using the pinning 
force at the center of the step F^ ter . Note that F^ ter is negative at low E ac increasing to 
zero at high E ac . We construct an adjusted pinning force F^l = F^ e + a x F^ ter . Note 
that Fpf^ = Fp^ e for symmetric steps. However, for asymmetric steps with a > 0, \F^\ is 
increased on the right edge and decreased on the left edge as compared to F^ e , yielding a 
more satisfactory fit to the data. The result is shown in Fig. 4 for the 2/1 and 1/1 steps 
where Sh ys t ~ \Fpf^\ x roughly describes the data. Here, instead of fitting the data to find 
the best values of a and A we have plotted the result for a = 5.0. 

Finally, regarding the discrepancies between our results and those obtained experimen- 
tally, we note that in all the simulations, for each given "sample", the largest hysteresis 
typically occurs when E ac = 0; in other words, mode-locking hysteresis occurred only in 
"switching samples." Furthermore, no chaos was observed for a wide range of parame- 
ter values. This is troublesome in the light of the experiments which show hysteresis in 
non-switching samples, and chaos at the step edges of switching samples. One possible 
resolution of this may come from further study of the depinning transition. Real CDWs 
exhibit two transitions, only one of which becomes hysteretic at low temperatures The 
Id model studied here exhibits one transition where a small hysteresis grows to a large one 
as a function of 71. However, there is evidence |16| that a 2d generalization of the model 
exhibits depinning transitions which more closely resemble the experiments. In that case 
a new "creep" motion arises between the first and second transitions, and perhaps this 
qualitatively different behavior will bare out these shortcomings. 
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In summary, we have studied the response of the global coupling model to a joint ac+dc 
field. Unlike the FLR model, this model exhibits hysteresis and nondegenerate mode-locked 
attractors. The model produces hysteresis on mode-locked steps exhibiting size asymmetry 
in the low ac-field regime. The pinning force and the step offset are the primary variables 
which control mode-locking hysteresis within a given phase diagram. Increasing rigidity of 
the CDW at large values of global coupling leads to low dimensional behavior and thus van- 
ishing of subharmonic steps. While the model exhibits no chaos in the regime studied, this 
low dimensional behavior may have some bearing on the low-dimensional chaotic attractors 
observed experimentally ||. 
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FIGURES 

FIG. 1. Mode-locking I-V curves for different values of the global coupling parameter, 71. Note 
that (a) is the FLR limit. The curves are shifted vertically for clarity. Some steps are labeled by 
their p/q values in (a). The emergence of asymmetric hysteresis with increasing global coupling is 
evident. 

FIG. 2. Mode-locking "phase diagram"; I-V curves for different values of E ac . Here 
k = 1.0, N = 256,7i = 2.5 and uj ex = 3.7. Curves are shifted vertically for clarity. The low- 
est curve is the regular I-V curve. Dashed lines connect the unlocking transition points for each 
step. Hysteresis varies in size for each step and is largest in the middle of the first oscillation. 
Note that the steps are offset to the right for low values of E ac and that the offset decreases as E ac 
increases. 

FIG. 3. Geometric representation of J c dw = E& c + F p in on an asymmetric mode-locked step. 
The step is drawn schematically to label important quantities and illustrate the nature of the 
asymmetry. 

FIG. 4. Hysteresis size vs. adjusted pinning force for 1/1 and 2/1 steps in the strong pinning 
regime. S hyst ~ |i^f + a x F^ ter \ x with a = 5.0 (see text). 
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